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ABSTRACT 

Objective: To investigate visceral leishmaniasis (VL) 
deaths and risk factors in two VL endemic areas of 
Bangladesh. 

Design: Retrospective cross-sectional. 
Setting: Two geographically and culturally different VL 
endemic subdistricts, Godagari in the district of 
Rajshahi and Trishal in the district of Mymensingh in 
Bangladesh, August 2009-December 2011. 
Participants: 51 094 inhabitants from randomly 
selected Unions in the two subdistricts. 
Main outcome measures: VL deaths, confirmed 
independently by qualified physicians using the verbal 
autopsy procedure ICD10 guideline. 
Results: The total number of people screened for VL 
deaths was 51 094 from 12 032 households from 
Godagari and Trishal subdistricts . About 16% of the 
people from Godagari were Tribals. The average age of 
the study population was 25.6 years (SD 18.4) and 
49.7% were females. The VL case fatality rate averaged 
6.12% (12/196) including 2/137 in Trishal and 10/59 
in Godagari. Most of the VL deaths (9/12, 75%) 
occurred at home and the rest in tertiary hospitals. 
None of these deaths had been reported in the national 
VL surveillance system. The VL case fatality rate in the 
Tribal ethnic (22.2%) population was about 17 times 
higher than that in the Bangali ethnic (1.3%) 
population (p<0.0001). Tribal ethnicity had an 18 
times (OR=18.1, 95% CI 3.6 to 90.6) higher risk for 
VL deaths compared with Bangali ethnicity (p<0.0001). 
Conclusion: VL deaths were found to be high in the 
study areas and were under-reported. The Tribal ethnic 
population was at the highest risk for VL deaths. The 
national VL Elimination Programme should give special 
attention to the tribal community in the endemic areas, 
especially for those in Rajshahi, and should strengthen 
VL surveillance by including tertiary hospitals in the 
national surveillance system. 
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INTRODUCTION 

Leishmaniasis is one of the most important 
vector-borne diseases with an estimated 1.98 
million disability adjusted life years and 
20 000-40 000 deaths per annum. 1 2 Recent 
empirical estimates found the annual 



Strength and limitations of this study 



This is the first study in the South East Asia 
Region, in which visceral leishmaniasis 
(VL)-related mortality was investigated in 
Bangladesh using the verbal autopsy procedure. 
The study is limited to only two unions of VL 
endemic districts. 

The recall bias in the study cannot be excluded, 
but we tried to minimise this by probing with dif- 
ferent events that occurred during the past 
2 years from the date of the interview. 



incidence of visceral leishmaniasis (VL) and 
cutaneous leishmaniasis (CL) to be 0.2-0.4 
and 0.7-1.2 million cases, respectively. 1 Ninety 
per cent of the global burden of VL occurred 
in India, Bangladesh, Sudan, South Sudan, 
Brazil and Ethiopia. 1 In the Indian subcontin- 
ent, the estimated population at risk for VL, 
annual VL incidence and VL death are 200 
million , 25 000 to 40 000 cases and 200 to 300 
deaths respectively. 3 For many years, there 
have been a substantial number of VL cases 
that have never been reported due to the lack 
of proper surveillance systems which contri- 
butes to the neglect of this disease. 

In Bangladesh, VL re-emerged in the early 
1980s, and since then more than 90 000 
cases have been documented. 4 5 Since 1996, 
the reported VL incidence in Bangladesh 
has ranged from 7000 to 9000 cases per year, 
but this is likely to be underestimated. 6 After 
adjusting the substantial under-reporting 
factor, the recent estimated VL incidence in 
Bangladesh ranges from 12 400 to 24 900 
cases per year. 1 The disease is prevalent in 45 
districts of Bangladesh; however, most of the 
cases are reported from the Mymensingh dis- 
trict. 4 5 The female Phlebotomas argentipes 
sandfly is the only vector; humans are the 
only reservoir and Leishmania donavani is the 
only species responsible for VL in the Indian 
subcontinent, including Bangladesh. 3-5 The 
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disease is fatal if not treated. Even with treatment, the 
mortality rate is high (10%).' In 2005, a national VL 
Elimination Programme has been launched jointiy by 
the Health Ministries from Bangladesh, India and Nepal 
to eliminate VL by 2015. 8 In Bangladesh, the goal of the 
elimination programme is to reduce VL cases below 1 
per 10 000 people in the VL endemic subdistricts 
(known as Upazila). 8 The strategies for VL elimination 
include active detection of VL and post-kala-azar dermal 
leishmaniasis (PKDL) cases, adequate treatment of VL 
and PKDL cases, interruption of disease transmission 
through integrated vector management, and social 
mobilisation. 8 9 Since 2006, several operational research 
activities have been conducted by government and non- 
government organisations to support the national elim- 
ination programme in Bangladesh. These programmes 
revealed that the passive case detection (existing surveil- 
lance system) highly underestimated VL cases in 
Bangladesh. 10-12 In 2008, a study conducted in the 
Rajshahi district showed that the incidence of VL was 27 
times higher than the elimination target. It also found 
that the knowledge of the community people about VL 
and its vector was poor and their health seeking behav- 
iour was also unsatisfactory. 12 

Similar to VL case reporting, it can be assumed that the 
mortality rate of VL is also highly underestimated in 
Bangladesh. Until now, there is no well-designed study to 
measure VL mortality in highly VL endemic areas of 
Bangladesh. According to the Directorate General of 
Health Services, the reported annual number of VL 
deaths varied from 6 to 36 during 1999-2009. 13 This 
appears to be highly underestimated because a 
community-based study in a small endemic area reported 
a high VL mortality rate (about 10.0%) in 2003. 14 The 
possible reason for the high underestimation of VL 
deaths might be a shortfall in the passive reporting 
system which does not include tertiary hospitals under 
VL case surveillance. VL patients who are poor and living 
in remote communities far from healthcare facilities 
often die at home without seeking treatment or prefer 
traditional healers for their healthcare and do not report 
to the public health facilities. Underestimation of deaths 
contributes to VL remaining a major neglected tropical 
disease. Thus, there is a need to investigate the magni- 
tude of VL deaths and its risk factors through a well- 
designed study and by using adequate tools for investiga- 
tion of death. Culturally, autopsy is not well accepted in 
Bangladesh. Verbal autopsy (VA) could be the alternative 
since it is the best available approach for assessing causes 
of death in the communities where most deaths occur at 
home. 7 Until now, no study has been conducted to inves- 
tigate VL deaths in the Indian subcontinent including 
Bangladesh after the start of the VL Elimination 
Programme. We have therefore undertaken this study to 
fill this knowledge gap in Bangladesh. The objective of 
this study is to estimate the VL mortality and risk factors 
in two geographically and culturally different and highly 
VL endemic areas of Bangladesh. 



METHODS 

Study sites and population 

The study was carried out in the two different VL 
endemic areas of Bangladesh, the Godagari subdistrict 
in the district of Rajshahi and the Trishal subdistrict in 
the Mymensingh district. Both Trishal and Godagari are 
highly endemic for VL. The study activities were per- 
formed from July to December 2011. Mymensingh dis- 
trict reports more than 60% of the total VL cases in the 
country. Trishal is one of the five most VL endemic sub- 
districts in Mymensingh and is located at 24.5792742°N 
90.3950388°E. Trishal consists of 12 unions and 160 vil- 
lages and all the unions are endemic for VL. Godagari is 
one of the highest VL endemic subdistricts in Rajshahi, 
which is located at 24.4667°N 88.3306°E with 9 unions 
and 396 villages and has a significant tribal population. 
Godagari was included in the study to investigate VL 
deaths among the tribal population and to compare 
them with those of the Bangali ethnic population. Based 
on our past experiences, we assumed that the tribal 
population was at higher risk for VL deaths. 

Sample size estimation 

It was reported that the VL case fatality rate (CFR) was 
about 10% in 2003 in VL endemic areas in Mymensingh 
district. 14 We assumed a significant reduction of CFR 
along with the reduction of VL incidence as a benefit of 
the VL Elimination Programme in Bangladesh. We 
therefore calculated the sample size assuming an 
average of 6.5% CFR with ±3.0% precision and 95% CI. 
A total of 191 VL cases were needed to obtain a 6.5% 
CFR in the study areas. To obtain a minimum of 191 
past VL cases in the past 2 years, the survey required a 
total of 47 750 people (considering an annual VL inci- 
dence of 20/10 000). The study areas included a popula- 
tion of 51 094 (24 934 from Trishal and 26 160 from 
Godagari). 

Study design and sampling 

The study design was retrospective cross-sectional. As 
mentioned above, we selected Trishal and Godagari sub- 
districts and one endemic union from each subdistrict 
was taken randomly (figure 1). The trained field 
research assistants (FRAs) visited each household, made 
household rosters and collected sociodemographic 
information from the head of the household. During 
their visit, FRAs also identified households with past 
cases of VL, any death and suspected VL deaths within 
the past 2 years. To minimise recall bias, the FRAs used 
probing with different events (mainly religious festivals) 
that occurred during the past 2 years from the date of 
the interview. We defined a suspected VL death as a 
death with a history of VL treatment or a death with an 
illness associated with two or more weeks of fever, plus at 
least one of the following: weight loss, enlarged 
abdomen, abdominal pain and/ or skin darkening. Then 
the VA method was applied to confirm the suspected VL 
death ( please see below) . A household with a VL death 
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Figure 1 Study profile (HH, 
household; VA, verbal autopsy; 
VL, visceral leishmaniasis). 



VL endemic areas in Bangladesh 



Endemic area with tribal 
community 



Endemic area without tribal 
community 



r \ 

Site-1 
Godagari Sub-district, 

Rajshahi District 
(5/9 endemic Unions) 



Randomly one Union 



Site-2 
Trishal Sub-district, 
Mymensingh District 
(12/12 endemic Unions) 

Randomly one Union 



Rishikul Union 
(6566 HHs from 55 Villages) 



Sakhua Union 
(5466 HHs from 7 Villages) 



HH screening and identification of the HHs with past history of VL and 
suspected VL death within last 2 years by Field Research Assistant 

— 196 VL patients (59 from Rishikul & 137 from Sakhua) 

— 174 VL HHs (51from Rishikul & 123 from Sakhua) 

— 382 deaths (174 from Rishikul & 208 from Sakhua) 

— 25 suspected VL deaths (22 from Rishikul & 3 from Sakhua) 



IfHH with past history of VL 



IfHH with suspected VL death 



HH head interview on socio- 
demographic status, knowledge 
about VL and health seeking behavior 
— 165 HHs (9 of 174 HHs were not 
interviewed due to migration out) 



VA by Physician for confirming 
VL death 

— 12 confirmed VL deaths from 
11 HHs (10 from Rishikul & 2 
from Sakhua) 



VL cured HHs (past h/o VL but 
not VL death in last 2 years) 
— 154 HHs (59 from Rishikul & 
137 from Sakhua) 



VL death HHs (past h/o VL and 
VL death confirmed by VA) 
— 11 HHs (9 from Rishikul & 2 
from Sakhua) 



Estimate VL mortality rate & its comparison between two sites, and 
identification of independent risk factors of having VL death 



confirmed by VA was defined as a VL death household. 
The household with a history of VL but no deaths in the 
past 2 years was defined as a VL cured household 
(figure 1). Household heads from both VL death and 
VL cured households were interviewed by the FRAs with 
a structured questionnaire for collection of information 
about the households' socioeconomic status, knowledge 
about VL and its vector, and for information about the 
treatment seeking behaviour used. At the same time, the 
coordinates (longitude and latitude) of all the house- 
holds in the study area were recorded using eTrex 
Venture Garmin single handheld GPS receivers. 

VA at suspected VL death household 

We developed an interview questionnaire for VA fol- 
lowed by an ICD 10 VA guideline and the questionnaire 
was reviewed by a VL expert. We also translated the ques- 
tionnaire into the local language, and it has been 
further reviewed, revised and approved by the icddr,b 
Ethical Review Committee. Briefly, the tool assessed the 



identity of the deceased and established the sequences 
of the events leading to death including signs and symp- 
toms of the illness before death. It is recommended that 
VA should be performed within 3 months to 2 years to 
minimise recall bias. 10 After identifying the household 
with a suspected VL death, a VL expert physician visited 
the household immediately with an FRA and conducted 
the VA. Respondents for the interview were household 
heads/ close relatives. Studies with VA showed that close 
relatives rarely refuse to give information and in some 
situations the response rates were 100%. 16-18 A second 
VL expert physician reviewed each VA questionnaire 
independently. In case of disagreement between the first 
and the second physician, the VA questionnaires were 
reviewed by a third expert physician and his verdict was 
final. 

Statistical analysis 

All data were entered and stored in a well-checked data 
entry program developed by EPI Info software (Centers 
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for Disease Control and Prevention, Atlanta, Georgia, 
USA). Data were cleaned and checked for duplicates. 
We explored the descriptive statistics to consider the 
nature of data. We calculated the CFR and compared 
the rates between different study sites (Godagari vs 
Trishal) and ethnic groups (Tribal vs Bangali) with the 
Z-test. We dichotomised socioeconomic, knowledge 
about VL, and health seeking behaviour variables for 
bivariate analysis. Variables which showed significant 
(p<0.05) association with VL death in the bivariate ana- 
lysis (Pearson's %2 and Fisher exact tests, where applic- 
able) were included in the logistic regression model. 
The final model was obtained by a backward selection 
strategy recommended by Kleinbaum et al. 19 The statis- 
tical software SPSS V.13.0 for windows (SPSS, Chicago, 
Illinois, USA) was used to carry out all the analyses. We 
constructed the study area map with spatial distribution 
of households in the study areas by using ArcView GIS 
V.3.2 (ArcGIS, ESPJ, Redlands, California, USA). 

RESULTS 

Sociodemographic characteristics of the study population 

A total of 51 094 people (12 032 households) were 
studied where 26 160 (6566 households) came from the 
Rishikul Union in the subdistrict of Godagari in the dis- 
trict of Rajshahi and the rest of the people from the 
Sakhua Union in the subdistrict of Trishal in the 
Mymensingh district. The average age of the study popu- 
lation was 25.6 years (SD 18.4) and 49.7% were females 
(table 1). About 16% (n=4243) of the population was 
tribal in Godagari, whereas all the people are non-tribal 
(Bangali ethnicity) in the Trishal area (table 1). 

VL-associated mortality based on VA 

During the period August 2009-December 2011, VL 
affected 0.4% (196/51 094) of the people in the study 
areas. The proportion of the VL affected population was 



higher in Trishal (0.5%, 137/24 934) compared with 
Godagari (0.2%, 59/26 160). Trained FRAs suspected a 
total of 25 VL deaths and VL expert physicians confirmed 
VL deaths in 12/25 by VA. The VL deaths in two areas 
averaged 6.12% (12/196) and accounted for 3.14% (12/ 
382) of the total deaths in the study areas. The VL CFR 
differed significantly from Godagari than Trishal and was 
17% (10/59) and 1.5% (2/137), respectively (p<0.0001) 
(table 2). The VL CFR in the Tribal population (22.2%) 
was about 17 times higher than the non-tribal/Bangali 
ethnic (1.3%) population (p<0.0001). 

Among the 12 confirmed VL deaths, 50% (6/12) were 
children aged below 17 years and 75% (9/12) were 
male. Most of the deaths occurred in tribal ethnic popu- 
lations, which accounted for 83.3% (10/12) of total VL 
deaths. About 33.3% (4/12) of all confirmed deaths did 
not receive treatment during their illness that led to 
death. Tuberculosis as an associated illness was found in 
8.3% (1/12) of the deceased. The majority of the deaths 
occurred at home and the rest were in the tertiary hos- 
pital (table 3). While checking the existing VL surveil- 
lance data, none of these 12 VL deaths were found. The 
spatial distribution of households with VL and VL 
deaths during August 2009-December 2011 is shown in 
figure 2. 

Factors associated with VL death 

To identify factors associated with VL death households, 
we compared household-related information between 11 
VL death households and 154 VL cured households. 
Bivariate association showed that households with a VL 
death differ significantly between Tribal-ethnic and 
Bangali-ethnic groups (p<0.0001). The VL death house- 
holds were socioeconomically poorer than VL cured 
households as daily labour household head, low monthly 
household income, and precarious housing structures 
were significantly more in VL death households, 
(table 4). The level of knowledge about VL was, however, 



Table 1 Sociodemographic characteristics and history of VL and death in the study population 





Trishal, Mymensingh 


Godagari, Rajshahi 


Total 


Total population 


24 934 


26 160 


51 094 


Total household 


5466 


6566 


12 032 


Mean age (SD) 


24.53 (18.91) 


26.70 (17.92) 


25.65 (18.44) 


Age ratio (years) 








<17 


10 868 (43.6) 


9309 (35.6) 


20 177 (39.5) 


17-45 


10 091 (40.5) 


12 895 (49.3) 


22 986 (45.0) 


>45 


3975 (15.9) 


3956 (15.1) 


7931 (15.5) 


Male population 


12 958 (52.0) 


13 267 (50.7) 


26 225 (51 .3) 


Ethnicity 








Bangali 


24 934 (100.0) 


21 917 (83.8) 


46 851 (91.7) 


Tribal 


00 (0.0) 


4243 (16.2) 


4243 (8.3) 


Past h/o VL in the past 2 years 


137 (0.5) 


59 (0.2) 


196 (0.4) 


Total deaths reported by respondent in the past 2 years 


208 


174 


382 


Total suspected VL deaths in the past 2 years 


3 


22 


25 


Total VL deaths in the past 2 years confirmed by VA 


2 


10 


12 



VA, verbal autopsy; VL, visceral leishmaniasis. 



4 



Huda MM, era/. BMJ Open 2014;4:e005408. doi:1 0.1 136/bmjopen-201 4-005408 



Open Access 




r-. 


CO 




CM 


CM 


CM 


CM 




CM 




lo 


CM 


CO 








LO 


CO 


CM 




CO 










CM 






CO 


0"> 


CO 


CO 


o 


CM 


CM 
















<M 






CM 


If) 


T 


CD 


r-. 


CM 


o 


O) 




UO 


O) 






T — 






T 




05 


ID 


oo 


CO 


CO 




^J" 


CD 


CM 




T 




LO 


O) 


oo 


CM 




CO 


CM 


co 


O 




CO 


o 


o 




CO 




CM 


CM 


CM 


CO 






UO 



O r-. 

d co 

CM t- 



O O) 
CO CM 



O CO 



i-i- O CM 



^ CO 
|v. CO 



i- CM i — CO 



CM LO 

CM O O) 

CM O CO 

CM i- 



lo cn 

■* 1- LO 



CO N O 

^ t- CO 

CM CD t— 

■* t- CO 

CM CM 



CO CO 

O 

C) T - " T~ 



o 



CD 







CO 


co 










CM 


CM 






TO 




!= 


2 


CCS 




CO 


CD 
> 




o 



almost the same between the two groups (table 4). 
Although the general health seeking behaviour of VL 
death and VL cured households was the same (p=0.431), 
their VL health seeking behaviour differed significandy 
(p=0.013; table 4). About 92.2% (142/154) of VL cured 
households visited qualified doctors when members of 
the family had VL, whereas only 63.6% (7/11) of VL 
death households visited qualified doctors for the same. 

The logistic regression model showed that ethnicity 
was the only independent risk factor for VL death. The 
adjusted model showed that households of Tribal ethni- 
city were at 18 (OR=18.1, 95% CI 3.6 to 90.6) times 
higher risk of having a VL death than households of 
Bangali ethnicity (p<0.0001) (table 5). 

DISCUSSION 

Our study is the first to estimate the actual mortality of 
VL in Bangladesh through VA, and the major findings 
of the study are as follows : VL mortality in the endemic 
community remains high, the tribal population is at 
higher risk for VL deaths and there is under-reporting 
of VL deaths by the existing VL surveillance system. The 
findings of our study are very important for the National 
VL Elimination Programme, because it may help the 
programme to make better plan to further reduce VL 
related mortality. 

Comparison with other studies 

It is alarming that the VL CFR is very high in the Tribal 
community in Rajshahi; it was found to be about 22%. 
The probable cause of the high VL CFR in that commu- 
nity could be their poor health seeking behaviours as 
they rarely visit government hospitals due to cultural 
beliefs and prefer traditional healers. 20 As a result, VL 
patients in that community remain untreated and die at 
home. We found that about 33% confirmed that VL 
deaths were not treated before dying and that all those 
untreated belonged to the Tribal community in Rajshahi. 
A previous study 12 conducted on the same community 
also identified a lack of access to proper healthcare, 
including high delay in VL diagnosis and treatment. 
Although Mymensingh is highly endemic for VL, the 
CFR of VL patients living in Mymensingh (all are of 
Bangali ethnicity) was not as high as that of those in 
Rajshahi district. Additionally, we found low VL mortality 
in Mymensingh as compared with the reported VL mor- 
tality in 2003. 14 The lower VL mortality in Mymensingh 
can be explained by the more focused attention paid to 
this highly endemic area by the national VL Elimination 
Programme of the country. 

The hospital-based VL CFR is substantially underesti- 
mated in the Indian subcontinent, including 
Bangladesh. 1 In Bangladesh, the existing hospital based 
surveillance system is mainly based on the passive report- 
ing system. There is no active follow-up of VL patients 
who had been treated in the hospital. As a result, the 
patients who died in the community/home are missed 
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Table 3 Characteristics of VL deceased confirmed by verbal autopsy 




VL deceased; n (%) 





Male 


Female 


uveraii 




M— a 


M— 1 
IN=o 




Age group (years) 








<17 


4 (44.4) 


2 (66.7) 


6 (50.0) 


17-45 


5 (55.6) 


1 (33.3) 


6 (50.0) 


>45 


0 (0.0) 


0 (0.0) 


0 (0.0) 


Ethnicity 








Tribal 


8 (88.9) 


2 (66.7) 


10 (83.3) 


Bangali 


1 (11.1) 


1 (33.3) 


2 (16.7) 


Received treatment before death 


5 (55.6) 


3 (100.0) 


8 (66.7) 


Having TB co-infection 


0 (0.0) 


1 (33.3) 


1 (8.3) 


Place of death 








Home 


7 (77.8) 


2 (66.7) 


9 (75.0) 


Hospital (tertiary) 


2 (22.2) 


1 (33.3) 


3 (25.0) 



TB, tuberculosis; VL, visceral leishmaniasis. 



by the existing surveillance systems. Until recently, the 
tertiary hospitals did not report VL cases to the national 
programme. Even very recently, only two tertiary hospi- 
tals started reporting VL cases to the national surveil- 
lance system. Our study results complied with these 
facts. Forty-two per cent (5/12) of VL deaths occurred 



at home, although they had been treated in the hospi- 
tals. Twenty-five per cent (3/12) died in the tertiary hos- 
pitals and the rest 33% (4/12) died at home without 
treatment. It is surprising that none of these 12 VL 
deaths were reported by the VL national programme 
even when the deaths occurred in the tertiary hospital. 



Rajshahi District 




Bangladesh 



Mymensingh District 




Mymensingh District 




Rishikul Union in Godagari Upazila, 
Rajshahi 




Legend 



• Tribal HHs with VL death in last 2 years 
a Bangali HHs with VL in last 2 years 
o Tribal HHs with VL in last 2 years 

Bangali HHs without VL in last 2 years 
Tribal HHs without VL in last 2 years 
I | Union boundary 



4 Kilometers 



N 

4- 



Sakhua Union in Trishal Upazila, 
Mymensingh 



^8, 



Legend 

• Bangali HHs with VL death in last 2 years 
o Bangali HHs with VL in last 2 years 
Bangali HHs without VL in last 2 years 
I Union boundary 



2 Kilometers 



N 

4- 



Figure 2 Study area map including spatial distribution of households, households with visceral leishmaniasis (VL) and VL 
deaths during August 2009-December 201 1 (HH, household). 
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Table 4 VL death and VL cured households and the bivariate analysis for socioeconomic, VL knowledge and health seeking 
behaviour-related information 





VI Heath HH 
vl ucctiii nn 


VI r-nroH HH 
vl luicu nn 






N-1 1 

IN— 1 1 


N-1 54 

IN— 1 iJt 






n \/o) 


n [/o) 


p ValUc 


Study site 








Mymensingh 


2 (18.2) 


115 (74.7) 


<0.0001 


Rajshahi 


9 (81 .8) 


39 (25.3) 




Socioeconomic status 








Illiterate HH head 


6 (54.5) 


74 (48.1) 


0.677 


Daily labour HH head 


10 (90.9) 


64 (41 .6) 


0.003 


Income about <US$ 100 


11 (100.0) 


101 (65.6) 


0.017 


Precarious house 


10 (90.9) 


83 (53.9) 


0.024 


Living room with <2 


7 (63.6) 


67 (43.5) 


0.223 


Family size with <5 members 


2 (18.2) 


57 (37.0) 


0.331 


Having bed-nets 


11 (100.0) 


152 (98.7) 


1.00 


Having domestic animals 


10 (90.9) 


130 (84.4) 


1.00 


Knowledge about VL 








Symptoms 


11 (100.0) 


136 (88.3) 


0.611 


Transmitted by mosquito/sandfly 


6 (54.5) 


71 (46.1) 


0.588 


The disease is curable 


11 (100.0) 


147 (95.5) 


1.000 


Health seeking behaviours 








Visit qualified doctor when get any sick 


1 (9.1) 


7 (4.5) 


0.431 


Visit qualified doctor when get VL 


7 (63.6) 


142 (92.2) 


0.013 


Ethnicity 








Bangali 


2 (18.2) 


127 (82.5) 


<0.0001 


Tribal 


9 (81.8) 


27 (17.5) 





HH, household; VL, visceral leishmaniasis. 



Thus, VL deaths are highly under-reported. Therefore, 
there is an urgent need for improvement of the existing 
surveillance system, implementation of active surveil- 
lance of VL patients under treatment and inclusion of 
the tertiary hospitals (especially those in the VL 
endemic areas) under a routine VL case reporting 
mechanism. 

Another important finding of the study is that about 8% 
of the VL deceased had tuberculosis, which indicates the 
existence of VL-tuberculosis coinfection in VL endemic 
areas of Bangladesh. The evidence of VL associated with 
tuberculosis has been reported in several VL endemic 
countries, including India. 21 However, this type of data has 
not been previously reported from Bangladesh. 

Although there are previous reports 22-24 on clinical 
risk factors for VL death, so far our study is the first 
where household-level risk factors were explored. 
Households with VL deaths were found to be 



socioeconomically poorer than those without VL deaths. 
Moreover, poorer VL health seeking behaviour was 
observed in households with VL deaths. Tribal ethnicity 
was the only independent risk factor for VL death in the 
adjusted regression model. This is because the tribal 
ethnic populations are comparatively poorer, less edu- 
cated and less aware about healthcare compared with 
the Bangali ethnic population. The VL Elimination 
Programme should therefore take the necessary steps to 
improve the VL health seeking behaviour of the tribal 
population, especially those in the Rajshahi VL endemic 
areas, to minimise their risk from VL death. 



CONCLUSION 

This study highlighted VL mortality in two VL endemic 
areas of Bangladesh, underlined the under-reporting of 
VL deaths and found tribal ethnic populations to be at 



Table 5 Independent risk factor on visceral leishmaniasis death (OR adjust using stepwise logistic regression) 


Variable VL death HH 


VL cured HH 


OR 


(95% CI) 


Adjusted OR* 


(95% CI) 


p Value 


Ethnicity 














Tribal 9 


27 


21.2 


(4.3 to 103.5) 


18.1 


(3.6 to 90.6) 


<0.0001 


Bangali 2 


127 


1.00 




1.00 






VL health seeking behaviour 














Unqualified doctor 4 


12 


6.8 


(1.7 to 26.4) 


4.3 


(1 .00 to 20.5) 


0.063 


Qualified doctor 7 


142 


1.00 




1.00 






'Adjusted by study site, household head occupation, income, housing structure. 








VL, visceral leishmaniasis; HH, household. 
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higher risk for VL death. A substantial number of VL 
patients never visit hospital for proper healthcare and 
die at home. The national VL Elimination Programme 
should give special attention to the tribal community in 
the VL endemic areas of Bangladesh and the existing 
surveillance systems must be improved by integrating 
active case detection as well as the active follow-up of the 
treated patients. 
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